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Using a Hamiltonian treatment, charged thin shells, static and dynamic, in spherically symmetric 
spacetimes, containing black holes or other specific type of solutions, in d dimensional Lovelock- 
Maxwell theory are studied. The free coefficients that appear in the Lovelock theory are chosen to 
obtain a sensible theory, with a negative cosmological constant appearing naturally. Using an ADM 
description, one then finds the Hamiltonian for the charged shell system. Variation of the Hamilto- 
nian with respect to the canonical coordinates and conjugate momenta, and the relevant Lagrange 
multipliers, yields the dynamic and constraint equations. The vacuum solutions of these equations 
yield a division of the theory into two branches, namely d — 2k — 1 > (which includes general 
relativity, Born-Infeld type theories, and other generic gravities) and d — 2k — 1 = (which includes 
Chern-Simons type theories), where k is the parameter giving the highest power of the curvature 
in the Lagrangian. There appears an additional parameter \ — (— l) fc+1 , which gives the character 
of the vacuum solutions. For X ~ \ the solutions, being of the type found in general relativity, 
have a black hole character. For x — ~ 1 the solutions, being of a new type not found in general 
relativity, have a totally naked singularity character. Since there is a negative cosmological constant, 
the spacetimes are asymptotically anti-de Sitter (AdS), and AdS when empty. The integration from 
the interior to the exterior vacuum regions through the thin shell takes care of a smooth junction, 
showing the power of the method. The subsequent analysis is divided into two cases: static charged 
thin shell configurations, and gravitationally collapsing charged dust shells (expanding shells are the 
time reversal of the collapsing shells). In the collapsing case, into an initially non singular spacetime 
with generic character or an empty interior, it is proved that the cosmic censorship is definitely up- 
held. Physical implications of the dynamics of such shells in a large extra dimension world scenario 
are also drawn. One concludes that, if such a large extra dimension scenario is correct, one can 
extract enough information from the outcome of those collisions as to know, not only the actual 
dimension of spacetime, but also which particular Lovelock gravity, general relativity or any other, 
is the correct one at these scales, in brief, to know d and k. 
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I. INTRODUCTION 

In a world with extra large space dimensions, of the order of some microns or less, as postulated in the 
gravitational field sees all the dimensions, whereas the standard model fields are confined to the usual world, or brane, 
of three space dimensions. Thus, the four dimensional spacetime of the brane can be seen as embedded in the d 
dimensional spacetime of the whole world. In this setup several different problems are tackled. For instance, one 
solves, in a way, the hierarchy problem, since both the Planck scale and the electroweak scale can be made of the 
same order, and so quantum gravity, being a electroweak scale phenomenon, can be experimentally tested. Indeed, by 
smashing particles against each other in experiments with the new generation accelerators or in eventual cosmic ray 
collisions, one can produce black holes, or perhaps other spacetimes with different causal structure, with tiny masses 
and radii. As this is a gravitational phenomenon, and the gravitational field spreads into all dimensions, the newly 
created spacetimes probe all the spacetime dimensions, and in addition render visible some quantum effects (see, e.g., 

Now, since in this setting one might expose quantum gravity and large extra dimension phenomena, in order to 
study the higher dimensional black holes, or other spacetimes with different causal structure, formed in the collision of 
particles, one should consider studying these objects in possible natural extensions to the theory of general relativity, 
conceivably appropriate to a further quantum framework development. The action of general relativity in d spacetime 
dimensions is proportional to the integral of the cosmological constant plus the Ricci scalar, and possible extensions 
should include powers of the Riemann tensor, such as the Kretschmann scalar, powers of the Ricci tensor, and powers 
of the Ricci scalar In principle one has to choose a criterion to pick up the right combination of curvature 
scalars. For instance, adopting the criterion of keeping the same degrees of freedom, one may wonder what is the 
most natural generalization of general relativity for other dimensions. Such a generalization is given by the Lovelock 
action [3], where essentially one keeps the field equations second order. In more detail, the Lovelock theory looks 
for an action which when properly varied yields symmetric tensors that are functions of the metric and its first and 
second spacetime derivatives, and divergence free. In four dimensions the corresponding action is the Einstein-Hilbert 
action, proportional to the cosmological constant and the Ricci scalar. When varied, the Einstein-Hilbert action 
yields the Einstein tensor and the metric tensor, the only tensors that display the above properties and contribute 
to the equations of motion. In dimensions higher than four Lovelock found that only certain precise combinations 
of higher powers of the curvature scalars could enter in the action pj. The interpretation in physical, geometric and 
topological terms of these precise combinations was put forward by Teitelboim and Zanelli [§] (see also 0, [IB]). In 
this context, one can first argue that in zero spacetime dimensions, a spacetime point, there is a zero dimensional 
topological invariant which is a numerical constant, called the cosmological constant. When this term is integrated 
in one or in two dimensions to form an action, i.e., when the term is dimensionally continued to the next odd or 
even dimension, it contributes to the equations of motion, which in this case are somewhat trivial, giving either a 
zero metric or a zero cosmological constant with an indeterminate metric. Now, in two dimensions, there is also the 
corresponding topological invariant. This invariant, the Euler characteristic, is obtained by integrating the so called 
Gauss-Bonnet term. In two dimensions the Gauss-Bonnet term is the two dimensional Ricci scalar, which can thus 
be consistently considered the corresponding Euler density. This term, being topological, does not contribute to the 
equations of motion in two dimensions. However, when both these terms, the cosmological constant (the topological 
invariant in zero dimensions) and the Ricci scalar (the topological invariant in two dimensions) are integrated in three 
or in four dimensions, i.e., when they are dimensionally continued to the next odd or even dimensions, they contribute 
non-trivially to the equations of motion, and indeed in three and four dimensions, yield the standard general relativity. 
Now, given general relativity in four dimensions, one can add a generalized Gauss-Bonnet (i.e., a generalized Euler 
density) term, in four dimensions, first discovered by Lanczos [llj, which when integrated in the four dimensional 
action gives essentially the corresponding Euler characteristic, also a topological invariant not contributing to the 
classical field equations. In turn, in five and six dimensions, not only an action with the cosmological constant and 
the Ricci scalar, mentioned above, contribute to the equations of motion, but also one can dimensionally continue 
the previous generalized Gauss-Bonnet term in order to have a meaningful extended action, which gives the desired 
equations of motion of second order. Of course, in six dimensions there is a new generalized Gauss-Bonnet term, 
which when integrated gives the corresponding Euler characteristic, not contributing to the equations of motion. 
Then, repeating this process of adding a generalized Gauss-Bonnet term of the previous even dimension to the next 
two dimensions, odd and then even, one gets the Lovelock gravity of that specific dimension. Thus, the Lovelock theory 
can also be considered as a dimensional continuation of the Euler characteristics of lower dimensions @, @, The 
theory has, in addition to Newton's constant and the cosmological constant, new arbitrary dimensionful parameters, 
which should be chosen carefully. We have given a classical motivation for including certain special new terms in the 
gravitational action in d dimensions, but there are also quantum motivations to study the dimensional continuation of 
the Euler characteristics. Indeed, in a string theory context, in order to prevent the existence of ghosts in a low energy 
limit of the theory, one has to add to the d dimensional general relativity action all the previous lower dimensional 
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Euler characteristics [l^. [T3|. 

Suppose now the world has indeed extra dimensions, and that the whole bulk spacetime obeys Lovelock gravity 
rather than Einstein gravity. Thus, in this case, it is important to study generic properties of black holes in the Lovelock 
theory. Solutions of spherical black holes in d spacetime dimensions for Einstein-Hilbert action plus a generalized 
Gauss-Bonnet term, i.e., a specially truncated Lovelock theory, were found in fl4L H5. fl^. flTl. fl8| . Thermodynamics 
were often studied in these works. In [13], where black hole solutions were also found, no relevant Lovelock term was 
put to zero, rather a special choice was made for all the Lovelock coefficients, which depend only on the cosmological 
constant and the dimension of the spacetime. All non-topological terms are included in the action, and the theory 
can thus be called dimensionally continued gravity. This choice is one among infinite, but it leads to a natural 
outcome, where the action in even dimensions has a Born-Infeld form, i.e., it may be regarded as the gravitational 
analogue of the Born-Infeld electrodynamics, and in odd dimensions has a Chern-Simons form. Also in [13], Lovelock 
theory was coupled to Maxwell electromagnetism, and the corresponding black hole solutions with electric charge 
and cosmological constant were found. Now, a natural extension of this prescription was advanced in [23], where 
it was allowed that the topological term plus some non-topological terms, starting from the top, could be put to 
zero. For instance, in d — 10 dimensions in this setting one may work just with the general relativity term and the 
following generalized Gauss-Bonnet term, by putting to zero the other higher dimensional Euler densities. For the 
nonzero terms, the choice of the coefficients of the theory is the same as in [lSj . This extension has as a special case 
the dimensionally continued gravity of [13]. By coupling this extension to Maxwell electromagnetism, a general class 
of charged vacuum spacetimes with a negative cosmological constant, including of course the black holes found in 
[13], were exhibited [23]. These vacuum solutions yield a natural division of the Lovelock theory into two branches, 
namely d — 2k — 1 > (which includes general relativity, Born-Infeld type theories, and other generic gravities) and 
d — 2k — 1 = (which includes Chern-Simons type theories), where k is the parameter that gives the highest power of 
the curvature in the Lagrangian. In the solutions there appears an additional parameter x, with x = (— 1) ■ This 
implies that, in addition, the solutions can be subdivided into two families of different characters, one family X = 1, 
comprising solutions of the type found in general relativity, has a black hole character (meaning that for a correct 
choice of the parameters, such as mass and charge, the solution is a black hole solution, although, of course, for other 
choices of parameters it can be an extremal black hole or a naked singularity), and the other family x = — 1, being 
of a new type not found in general relativity, has a naked singularity character (meaning there is no possible choice 
of the parameters that gives a black hole solution, the full vacuum solution is always singular without horizons). Of 
course the solutions also include empty spacetimes. Since in general there is a negative cosmological constant the 
spacetimes are asymptotically anti de Sitter (AdS), or when empty they are AdS itself or some form of it. 

Given that one has a family of spacetimes with either black hole or naked singularity character, it is now important 
to study matter effects, either static or dynamic (collapsing or expanding), on these solutions. Such a phenomenon 
could be interesting in the aftermath of a collision between charged particles where the debris formed in the collision 
can be accreted or excreted in the newly formed charged spacetime. Accretion and excretion processes are usually 
technically elaborate, so to turn the analysis simpler one can study the case of a thin shell in the background spacetime, 
be it a black hole or otherwise. Static and dynamic solutions on such backgrounds are then of interest. However, 
in Lovelock theories, even a thin shell and its dynamics can bring complications. Indeed, in general relativity for 
instance, to study a thin shell in a g iven background one has to make a smooth junction from the interior to the 
exterior spacetime, as was done in [2JJ] . Now, the junction conditions, in general, depend on the theory one is studying, 
and for theories with higher powers on the Riemann and other tensors this can be nontrivial [13] ■ An elegant and 
useful approach, that bypasses in a way the junction conditions for a spacetime with a thin shell, uses a Hamiltonian 
formalism for the theory under study. This approach was developed by Hajlcek and Kijowski [23] for the theory 
of general relativity with matter in four dimensions, and was subsequently explored by Crisostomo and Olea in d 
dimensional general relativity coupled to Maxwell theory and matter with applications mainly in three spacetime 
dimensions [24 . l25| . The method requires that one puts the theory in a Hamiltonian form and then one directly 
integrates the canonical constraints, producing the shell dynamics for the background one wants to study. The pay 
off is that, not only one avoids to have to develop covariant junction conditions for each particular theory one is 
studying, but moreover the full spacetime, comprised of interior, shell, and exterior components, is treated as whole. 
This method is particularly well suited for symmetric configurations. The great power of the formalism was brought 
up in [26j where it was applied in pure Lovelock gravity to thin shell collapse in the background of the uncharged 
spacetimes found in [13] ■ Now, since the newly formed spacetime, black hole or otherwise, generated from a collision 
between charged particles, is most probably charged, and for the same reason the collapsing or expanding debris are 
also charged, it is of interest to study, using the Hamiltonian formalism advanced in [26j, the gravitational dynamics 
of a charged thin shell in the charged sector of the spacetime solutions, black hole or otherwise, found in [13], of 
Lovelock gravity coupled to Maxwell electromagnetism. In the usual extra dimension scenarios the electromagnetic 
and matter fields, being confined to the brane, do not probe the extra dimensions, so that an axisymmetric shell, 
static or collapsing, around a spherically symmetric black hole is an important configuration to study. However, in 
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order to simplify the analysis we study instead, as a zero order approximation, a spherically symmetric shell around 
a spherically symmetric black hole (such a configuration would be quite realistic for charged fields and particles that 
can probe the extra dimensions). Of course, within general relativity, a very special case of Lovelock gravity, such 
an analysis should recover the earlier results, obtained through a junction condition formalism [2l[ of charged shell 
dynamics in a charged black hole spacetime [27I. \2tL |29| . All the works just mentioned are classical, and an extension 
of full Lovelock gravity into the quantum domain, or even within a semiclassical approximation, although important, 
seems beyond reach. 

In the present article we extend the treatment of the classical dynamics of an uncharged thin shell in an uncharged 
spacetime, black hole or otherwise, with a negative cosmological constant background in Lovelock-Maxwell theory, 
given in [2y], to the classical dynamics of a charged thin shell in a charged spacetime, black hole or otherwise, with a 
negative cosmological constant background in Lovelock-Maxwell theory, and among other things, show the power of 
the Hamiltonian method in unifying in a single unit the three sectors of the problem, namely, the gravitational, the 
electrodynamic and the matter sectors, as well as taking into account automatically the smooth junction of the several 
regions of spacetime in question, the interior, the exterior and the thin shell in between. In Sec. [ill the system under 
study, namely charged vacuum plus charged thin shell, is initially put in an action and Lagrangian framework, in both 
tensor and differential forms languages. The field content of the action is composed of three sectors, the gravitational, 
the electrodynamic, and the matter term. In the gravitational action we establish the coupling constants and make the 
choice of Lovelock coefficients, cutting off the Lovelock polynomial at the highest possible power of the curvature for 
which the theory is sensible. In the electrodynamic action we write the Maxwell term and add a current term, which 
describes the current of charged matter in the shell, and precede each term with its coupling constant. In the matter 
action we define the energy-momentum tensor of a perfect fluid matter source. We then write the same action in the 
Hamiltonian form in tensorial language. For each of the above mentioned sectors of the action we define the canonical 
coordinates and its conjugate momenta, up to surface terms. As the system under study is a system with constraints, 
we write the constraints and their respective Lagrange multipliers for each of the sectors of the action. Next, we 
obtain the Hamiltonian field equations from the Hamiltonian action by using spherically symmetry in the action and 
varying it with respect to the canonical coordinates and momenta, for the evolution equations, and the Lagrange 
multipliers, for the constraint equations. In Sec. IIIII we analyze the vacuum solutions of the derived Hamiltonian 
equations and describe their properties. Afterwards we establish the geometrical framework of the thin shell and its 
Hamiltonian description. We divide spacetime into three parts, interior spacetime, thin shell, and exterior spacetime, 
and describe their geometric set-up. Then we specify the matter properties of the thin shell, which we define as being 
that of a perfect fluid. We write its energy-momentum tensor and its relevant projections. We then solve the complete 
Hamiltonian equations around the thin shell. First, we note that the vacuum solutions yield a natural division of the 
Lovelock theory into two branches, namely d— 2k — 1 > and d— 2k — 1 = 0, where k is the parameter that gives the 
highest power of the curvature in the Lagrangian. There appears an additional parameter x, with x = (~ l) fe+1 , which 
gives the character of the solutions, namely, the vacuum solutions may have a black hole character (x = 1) or a naked 
singularity character (x = —1). Since there is a negative cosmological constant, the spacetimes are asymptotically 
AdS, and AdS when empty, or some form of it. Second, the whole integration of the equations from the interior to 
the exterior vacuum regions through the thin shell is performed, where the smooth junction of the several regions of 
spacetime in question is automatically taken into account by the integration, showing definitely the efficacy of the 
method. The shell's dynamics in the vacuum spacetimes is then derived and the electrodynamic constraint equation 
recovers the electric charge conservation. An analysis of the equations of the thin shell is performed by studying two 
interesting cases, namely a static thin shell in equilibrium and the gravitational collapse of a thin shell (the study of 
gravitational expansion is simply the time reversal of gravitational collapse, and so it is not necessary to analyze it in 
any detail) . For the static shell we determine the radii at which the shell is in equilibrium and the pressure necessary 
to maintain the shell at this radius. For the gravitational collapse of a thin shell we start studying a simple example 
of dust matter in an empty interior and then prove cosmic censorship in a general case. More specifically, following 
the natural division of the Lovelock theory into the two branches mentioned above, we study the collapse of a thin 
shell onto an empty interior without cosmological constant, and give as examples collapse in four and ten dimensions 
in general relativity, collapse in ten dimensions in Born-Infeld type theories, and collapse in ten dimensions in general 
relativity with a single generalized Gauss-Bonnet term. We also study gravitational collapse in five dimensions onto an 
empty interior with a cosmological constant as an example of Chern-Simons type theories. Plots are drawn for these 
cases just mentioned, which illustrate the dynamics of the collapse and exhibit the cosmic censorship at work. This 
example of collapse of a shell into an empty interior, in a number of theories specified by d and k, shows that electric 
charge provides a mechanism for cosmic censorship for all relevant cases, including those which, in the uncharged 
case, did not conform to it. Then, more generally, for gravitational collapse onto black hole interiors, it is proven that 
the equations respect the cosmic censorship hypothesis. In Sec. [lV]we conclude and draw some physical implications 
in connection with the extra large dimensions scenario. We put c = 1. 
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II. ACTION, LAGRANGIAN, HAMILTONIAN AND EQUATIONS OF MOTION IN LOVELOCK 
GRAVITY COUPLED TO MAXWELL ELECTRODYNAMICS AND A CHARGED THIN SHELL IN A 

SPHERICALLY SYMMETRIC BACKGROUND 

A. Action, Lagrangian, and Hamiltonian in a general form 

1. Action and Lagrangian 

For our purposes of studying charged matter in a d dimensional charged background spacetime, the field content 
is divided into three sectors, the gravitational, the electrodynamic, and the matter sectors. Then, the action / can 
be written as the sum of a gravitational action I^ 9 \ an electrodynamic action I^ e > , and a generic matter field action 
/( m ), which will be specialized later on to be a charged thin shell. Thus, 

I = I {9) + I {e) + I (m) . (1) 

The gravitational action and Lagrangian: 

The gravitational sector of the action, depends on which theory one wants to adopt. General relativity, first 
formulated in four spacetime dimensions, can be trivially extended to higher d spacetime dimensions by changing the 
action from a four dimensional integral of the cosmological constant term and the Ricci scalar, to a d dimensional 
integral of both terms. However, this extension is no longer unique. Another natural extension is given by the 
Lovelock gravity Q, in which its action results from demanding that the Euler-Lagrange equations derived from the 
corresponding Lagrangian yield all tensors symmetric in fiis, occurring concurrently with the metric tensor 
and its first two derivatives, i. e., — ^4 MI/ (g p(J ; g p(T;7 ; g p(7 :ys), and divergence free. In four dimensions only the 
cosmological constant and the Ricci scalar yield an action with these properties, which is precisely the Einstein-Hilbert 
action. As one goes to higher dimensions new generalized Gauss-Bonnet terms, topological in nature in the previous 
dimension, make their appearance. As mentioned in the Introduction, the Lovelock theory can then be considered as 
a dimensional continuation of the topological Euler characteristics of lower dimensions @, H, The theory has, in 
addition to Newton's constant and the cosmological constant, new arbitrary dimensionful parameters, which should 
be chosen carefully. Due to its interest and generality, we work here with Lovelock gravity. The gravitational part 
1^ in (Q]) is then the Lovelock action, which turns out to be a polynomial in the curvature tensor, of degree [d/2], 
where the brackets [d/2] represent the integer part of d/2. This action, and the corresponding Lagrangian £ defined 
through I = J d d xC, can be written as 

fr) =k d d xJ2 «p a/=I2-> C-'-'C K% ■ ■ ■ RZ-\Z > ( 2 ) 

J M p = 

where k is inversely proportional to the Newton's constant (which will be appropriately chosen below), M. stands 
for the spacetime manifold, g is the determinant of the spacetime metric, and /ii,^2, ■■■ are spacetime indices. The 
generalized S function is antisymmetric in all of the upper indices and all of the lower indices, and i?p lP 2 2 is the 
Riemann tensor. The coefficients a p are arbitrary in general, apart from the first two ao and a\. Indeed, in the 
gravitational action (J2j) , the first term of the integrand is the cosmological constant A and the second term of the 
integrand is the Ricci scalar R, i.e., the terms of the Einstein-Hilbert action. This shows that general relativity is 
contained in the Lovelock theory as a particular case, namely by putting in the action all the a p — for p > 2. 
For even dimensions, the term p = d/2 in the action |f5]) is the Euler characteristic of that d dimensional manifold 
and does not contribute to the field equations. However, although not relevant for the purposes of this paper, the 
presence of the Euler term guarantees the existence of a well defined variational principle for asymptotically locally 
AdS spacetimes (see [IS])- Lovelock gravity, from the way it is constructed, has the same essence as general relativity. 
However, the theory is more complicated, it sets extra problems not present in general relativity and yields new 
interesting features. For instance, for d > 4 the fields may evolve in a non-unique manner, such that, given initial 
values for the fields at t = to, at t > to the equations of motion do not determine those fields completely. This 
is due to the presence in the Lagrangian of high powers in the first derivative of the metric tensor with respect to 
time [ala, llflll- An interesting aspect of Lovelock gravity is that, although its linearized approximation is classically 
equal to the corresponding linearized approximation in general relativity in the full strong gravity regime of the 
theory, the higher powers of curvature in the Lagrangian yield solutions that are different, and such type of solutions 
in Lovelock theory cannot be reached through a solution in general relativity 

Another feature is the fact that, in addition to the constant k (inversely proportional to an appropriate generalization 
of Newton's constant) and the cosmological constant A, there are the [(d+ l)/2] arbitrary dimensionful parameters. 
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Now, for a given dimension and an arbitrary choice of the coefficients a p , the dynamical evolution can become 
unpredictable [1, 0, [T3| , thus it is advantageous to restrict these coefficients in order to be able to construct meaningful 
black hole and other solutions. An important step in this direction was taken first in [T{| and generalized in [2(3], where 
the coefficients a p were chosen so that unique sensible solutions, with well defined perturbations could be found. One 
way to meet these requirements is to demand, in a consistent way, that the theory has a unique cosmological constant 
A as shown in 0, H3|j which leads to the choice 

{ ^(p-k) ( k\ , 
Ta=m\p)> P^ k , (3) 
0, p > k 

where I is a length scale of the theory given in terms of the cosmological constant A by 

A = -<t«. (4) 

As seen in the last expression, in this setting the cosmological constant is negative throughout, and so the solutions 
will be asymptotically AdS solutions. We are then left with a set of theories labeled by an integer k, with 1 < k < 
[(d — l)/2]. One should explain better the dependence of the theories on k. In [llj k was put equal to its maximum 
value always, k = [(d — l)/2], i.e., no pertinent Lovelock term is put to zero, and one has the corresponding special 
choice made in Eq. $3§. In this case, all non-topological terms, plus the topological non-dynamical term in the 
even dimensional case, are included in the action, and the theory can thus be called dimensionally continued gravity, 
since as one goes one dimension up, what was a topological term in the previous dimension, is now a term that 
contributes dynamically when one continues the dimension. In even dimensions, d > 4, this choice leads to a Born- 
Infeld gravitational action, the gravitational analogue of the Born-Infeld electrodynamics, and in odd dimensions, 
d > 5, it leads to a Chern-Simons gravitational action This dimensionally continued theor y fTgj ] can be naturally 
extended by canceling higher order terms, from a given integer k upwards. This was done in [20j. For instance one 
can set up a theory in which only the terms p — and p = 1 appear, thus k = 1. This theory is general relativity. If 
one puts k = 2 one gets general relativity plus the corresponding generalized Gauss-Bonnet term, and so on, up to 
the dimensionally continued gravity where k = [(d — l)/2]. In brief, for a given dimension d, depending on the choice 
of the integer k, one generates different theories, with k giving the highest power of the curvature in the Lagrangian. 
Now, the other fundamental constant k can be related to a generalized Newton's constant Gk, labeled by k, with 
units [Gk] = ( length) d ~ 2fc , through 

(0 



2{d-2)\n {d _ 2) G k ' 

where fld-2 is the area of the unit sphere in d— 2 dimensions. These theories have well defined black hole configurations. 
The electrodynamic action and Lagrangian: 

The electrodynamic sector of the action (JTJ) , the electrodynamic action I^ e \ can be chosen to be the electromagnetic 
Maxwell term plus a matter electric current, and can be written as 

/ (e) --^^/ d d a;V=iV u ' -- / d d x^J»A Pl (6) 



where e is related to the vacuum permitivity eo through the expression e = jff-^ , is the determinant of the spacetime 

metric Q^, F 2 = F^F^ is the square of F M „ = V^A^ — V^A^, which is the Maxwell tensor, with A p being the 
electromagnetic potential, and is the electric current. For vacuum solutions J M is zero, and one has pure Maxwell 
electromagnetism. 

The matter action and Lagrangian: 

The last term in the action |T]), the matter action /^ m ' is defined as 

(7) 

where the matter Lagrangian C m is defined through the energy-momentum tensor of the matter by 

(d-2)!» (d _ 2) SIM 

where 5 here means variation. In vacuum /( m ) = 0. Later on we will specify the matter as being made of a charged 
thin shell. 
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2. Hamiltonian 



In order to apply the Hamiltonian framework, we use the Arnowitt-Deser-Misner (ADM) formulation [30( where 
there is a foliation of spacetime into t = constant hypersurfaces, denoted by E t . In this foliation of the spacetime the 
metric, both inside and outside, is written generically as 

ds 2 = -{N^ fdt 2 + g t] {N l dt + dx l ){N ] dt + dx J ) , (9) 

where N 1 - and N l are the lapse and shift functions of the foliation, respectively. The <7y are the metric coefficients of 
the intrinsic geometry of the hypersurfaces E t , where i,j run only on the spatial components. It is known that the 
action I — J d d xC, with C the Lagrangian, can be written as 

1= J dt J d d ~ x x (n^glj - H) , (10) 

where ir t: > = SC/Sglj are the momentum components conjugate to the metric components gij, gij are the respective 
coordinate time derivatives, and H is the Hamiltonian of the system. It is then useful to write the Hamiltonian in the 
form 

H = N ± H±_ + N*7U , (11) 

where H± is the normal Hamiltonian constraint that generates the time translations normal to the hypersurface E t , 
and Hi are the tangential constraints that generate the translations in each of the hypersurfaces E t , which is the same 
as saying that is the generator of hypersurface diffeomorphisms, or that it generates coordinate transformations 
in Sj. In addition, one can write, 

n± = n { l ] + h { 1 ] + H { r } , "Hi = n\ g) + + n ( r ] (12) 

where the superscript (g) refers to the gravitational sector, the superscript (e) refers to the electrodynamic sector, 
and the superscript (to) refers to the matter sector of the Hamiltonian. The lapse function N 1 - and the shift functions 
N 1, of the ADM metric act here as Lagrange multipliers. 

The gravitational Hamiltonian: 

For the Lovelock theory, the Hamiltonian components of the gravitational field can be taken from the action j2]| , the 
momenta conjugate to the corresponding Lagrangian 7T 1 ^ = SC/Sglj, and glj, with the result [si. I9L [Io| 

n { [ ] = -KJgY (d ~ 2g)! " p Bit ■ ■ ■ B? p -T , (13) 

P 

= -2^.. (14) 

The Rijki in (fl3|l are the spatial components of the curvature tensor in the d— dimensional spacetime. The Gauss- 
Codazzi equations give us the relation between this tensor and Rijki, the Riemann tensor intrinsic to the surface 

Rijki — Rijki + K ik Kji ~ KuKjk , (15) 
where is the extrinsic curvature tensor of the surface, given by the expression, 

Kij - C^lj Wid + N Jti - gij) . (16) 

The semicolon in the indices denotes the intrinsic covariant derivative in the hypersurface Ef, and the dot is the deriva- 
tive of the spatial-spatial components of the metric with respect to the time coordinate. The conjugate momentum 
to g^ can be found through SC/Sglj, which in this case gives [I. fsL flo| 

p— 1 

i P- (d~ 2p)!a p n x iii-ia.-i2p-i „ 
*J = l; \ 7 ll^ 2^+1 2^»(ri*wW a ."-i£-i x 

p s— 

v . . . jihs-ite, K i2s+i Tftep-i / 17 \ 
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where, 



Here the double factorial is defined by: 



(_4)P-2 

° s ^= s![ s (p (18) 



s- (s-2)...5-3 - 1 s > odd 
s!! = <( s • (s - 2) . . .6 • 4 • 2 s > even . (19) 
1 s = -l,0 

The electrodynamic Hamiltonian: 

The action in a Hamiltonian form of the electrodynamic field in a curved background is given by I e = 
J dt J d D ~ 1 x p l Ai — ^N 1 - (ri d - 2 ^P t Pi + an^a ^'-^J") + ~ ^ J° > where p l is the momentum conjugate to 
the spatial components of the gauge field Ai, ip= Aq, py are the spatial components of the Maxwell tensor, J° is the 
time component of the electric current, and 11^—2 is the area of the (d — 2)-dimensional unit sphere. We are ignoring 
surface terms since they will be automatically taken into account in this setting. From the electrodynamic action, 
and using the definition given in Eqs. ifTTj) and (f!2| . we have 

V.f = \^ 2 ^p l + -^-F^F l0 , (20) 

Vlf = 0. (21) 
Besides the usual constraints, there is also a new constraint E v , associated with the Lagrangian multiplier ip 

E v = p\ t - J° . (22) 

The matter Hamiltonian: 

The matter sector of the Hamiltonian constraint is written as 

«i m) = V9Txx, (23) 
n\ m) = 2^/gT u , (24) 

where is the matter energy-momentum tensor, and where the index _L means that the tensor has been contracted 
with the hypersurface normal n M = (— iV- 1 , 0, 0, 0). In the following we will specify the matter as being made of a 
charged thin shell. 



B. Hamiltonian and field equations for thin shells with interior and exterior static vacua in spherically 

symmetric spacetimes 

1. Gravitational, electrodynamic, and thin shell Hamiltonians in spherically symmetric spacetimes 

Using the Hamiltonian formalism, we now want to set up the field equations appropriate for charged thin shells in 
static spherically symmetric charged Lovelock backgrounds. The formalism advanced here can be used for both static 
and dynamic thin shells. In the next section we will make full use of the formalism when we apply it first to find the 
pressure necessary to maintain a charged thin shell in static equilibrium in a vacuum background interior, black hole 
or otherwise, and second to the study of the collapse of a charged thin shell into a vacuum background interior, black 
hole or otherwise, in Lovelock gravity coupled to Maxwell electromagnetism. Since the three sectors that enter the 
problem are the gravitational background, the electrodynamic interaction, and the matter that constitutes the shell, 
we have to set up the metrics for the interior and exterior to the shell, we have to give the form of vector potential 
field, and we have also to give the energy momentum tensor. 

The gravitational Hamiltonian in spherically symmetric spacetimes: 

We assume now the shell traces its spacetime trajectory in a static spherically symmetric background. Thus the 
generic ansatz both for the exterior and interior to the shell can be written as 

dr 2 

ds 2 = -N 2 (r)f 2 (r)dt 2 + —— + r 2 dVL 2 d _ 2 , (25) 
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where Schwarzschild type coordinates {t, r, 9 1 , 9 2 , ■ ■ ■ , 9^ d 2 )} have been chosen, N(r) is the lapse function exclusively 
dependent on the radial coordinate, f(r) is the metric function, dependent on r only, and dil 2 _ 2 is the line element of 
the d - 2-dimensional unitary sphere, written explicitly as dVt 2 d _ 2 = (dd 1 ) 2 + sm(9 1 ) 2 (d9 2 ) 2 + s'm(8 1 )' 2 sm(9 2 ) 2 (dd 3 ) 2 + 
h nf=i sin(0') 2 (d9 < - d - 2 y) 2 . With the ansatz {25}, the first gravitational constraint (Taj) is written as 

and the other gravitational constraints {H| yield the equation 

,(<?) /om o^y* p!(d-2p)!a p 

2P+ 1 (d-2p- I,. 

p=Q y r ' s=0 



= K Vn {d - 2) ! g ^ _ g y^/x-rt (i - fy x 



x ^( rd " 2p " 1 K f ) 2p_2s " 1 ) ' ( 2? ) 

where is the angular part of the determinant g of the intrinsic metric gij of the hypersurface St, with g — g rr r d ~ 2 Q, 
and we use the fact that the metric is diagonal. The function D s /p\ is defined in Eq. fT8|) , 

The electrodynamic Hamiltonian in spherically symmetric spacetimes: 

For a static system with spherical symmetry the electromagnetic field is purely electric and radial for an outside 
observer at rest. In this case, the vector potential has only one non-zero component, which depends exclusively on 
the radial coordinate. 

A t = A(r) . (28) 

Given the symmetries and the fact that we are working on a static background, the totally spatial components of 
the Maxwell tensor are null, i.e., — 0, i,j = r, 6 l , ■ ■ ■ , 9^ 2 \ The only non-vanishing component of the Maxwell 
tensor is F tr — —d r At(r), with also F rt — d r A t (r). Thus, electrically charged, static, spherically symmetric vacuum 
solutions imply 

Fy = 0, (29) 
p* = (0,^,0,..., 0) , (30) 
Ai = and f = , (31) 

where Eq. lj29j) means there is no magnetic field, Eq. ([30"ll means the electric field is spherically symmetric, and Eq. 
([3T]l means the field is static. From equations (|20"ll - lf2"2"j) . the electrodynamic constraints, associated with the Lagrange 
multipliers N 1 - and ip, are, respectively, 

n [ l ] = ^=^-2PV = ^=^- 2 b r ) 2 r 2 M, (32) 

= 0, (33) 
E v = p r r - J° , (34) 

where /(r) comes from |25|) . 

The thin shell Hamiltonian in spherically symmetric spacetimes: 

We want to spell out completely the matter constraints, namely, — ^/gT±± and H- m ' ) = 2^fgT^\ for a charged 

thin shell in a charged vacuum background geometry. More precisely we want to develop a Hamiltonian framework 
for this situation and set up a natural junction which arises in this formalism. For this to be made we have first 
to describe the geometrical set-up in order to be able to write the projected components of the energy-momentum 
tensor, namely, T±± and Tj_j. 

First we state the nomenclature for the thin shell and for the interior and exterior spacetimes. The thin shell is 
ad — 2-dimensional spacelike surface, which evolves in time, see Figure [TJ This time evolution of the thin shell can 
be represented by a timelike hypersurface dV^, a boundary surface, which divides spacetime into two regions, the 
interior, denoted by V^', and the exterior, denoted by V^ + \ At each point on the boundary surface there exists 
a unit space-like vector £, with components £ M , normal to dV^, and pointing from the interior to the exterior 
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In dVz there is a set of intrinsic spacetime coordinates p a , where a runs from to d — 2. In the regions 
and y(+) independent coordinates are introduced, and respectively (where p runs from to d— 1), and so the 
parametric equations for dV^ in these coordinates are xt(p a ) and x^.(p a ). Thus, 9Ve has tangential vectors e a with 
components given by = Since represents the spacetime evolution of the d — 2-dimensional thin shell, the 
d- velocity u, with components u^, of the matter of the shell is tangential to dV^, with u being orthogonal to £, and 
vanishing outside the shell. One can also consider an intrinsic velocity vector u, with components u a , related to the 
velocity u through the relation u M = e%u a , where u is built by considering the evolution of matter in the shell using 
coordinates (p 1 , • • • , p d ~ 2 ) intrinsic to the shell, and the intrinsic time p°. 



PIG. 1: The shell trajectory manifold 9Vj in spacetime is depicted. V_ and V+ are the spacetime regions interior and exterior 
to <9Vj, respectively, and E t is the hypersurface intersecting dV%. Also shown are the normal vector n (normal to Et), the 
velocity vector u (tangent to dV$), the vector £ (the spacelike normal to dV^), and the vector e a (the generic tangent vector to 
dV{). The vectors indicating the origin and direction of the adapted coordinate system {T, X} are depicted together with the 
u and £ vectors. 

Working with the spherically symmetric metric (|25f , we choose the intrinsic coordinates of the thin shell surface 
dV^ as the proper time r, and the spherical angles of the Schwarzschilddike coordinates, so that p a = (r, 9 1 ,..., 6 d ~ 2 ). 
The motion of the shell is governed by the radial function of the proper time r — R(t). The derivative with respect 
to r will henceforth be denoted by a dot, e. g. 4zR{f) = R(r). The line element for dV^ in these coordinates is 



We want to study the spacetime evolution of a thin shell in a vacuum background with metrics of the type given in 



The radial and angular coordinates match continuously at the shell, but the time coordinates do not in general. The 
radial and angular coordinates match on account of the fact that the surface area of the shell is independent of the 
coordinates used to determine it, and it should be the same on both sides, as we are working in the limit of infinitesimal 
thickness, the thin shell limit. As the area of the shell is a function of the radial position of the shell R(t), there 
should not be a different radius for the shell inside and outside, as we would be attributing two different areas to the 
shell. Also, spherical symmetry allows us to induce coordinate charts on the manifold dV^ merely by restricting the 
charts of and to their boundary on dV^. As, due to symmetry, the radial and angular coordinates on both 
sides of the shell are the same, their restriction to the shell should be the same, and hence match perfectly. However, 
the restriction of the inside line element (|36|) to the surface of the shell must match the restriction of the outside line 
element f37|) . all having the general form f35|) on the shell. In these coordinates the vectors u and £ defined above 
(see also Figure [J) have the following components, 




dr 2 + R 2 (r)dn 2 d 



(35) 




(36) 
(37) 




(38) 
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C = 1-^,7,0,0) • (39) 



where 



where j± is a generalized Lorentz factor. These vectors were obtained through the matching of |36|) and l(37|) to the 
surface line element ([35| , using the relation u M = , where the t± coordinates and the radial coordinate are treated 



7=V/ 2 + ^ 2 > ( 4 0) 



dr 

as functions of the proper time r, We have also used u^u^ — —1, £ M it^ = 0, and = 1 (as £ is the spacelike normal 
todVt:). 

The shell matter properties are described by the surface energy- momentum tensor T ab , which is orthogonal to £ 
and vanishes outside the hypersurface dV^. For an observer in the shell, T ab is written in intrinsic coordinates, which 
for a perfect fluid is given by 

T ab = <7u a u b + P (h ab + u a u b ) , (41) 

where a is the rest mass surface density, P is the surface pressure (or, when negative, the surface tension), and 
h ab = g^e^e^. Tab is confined to the hypersurface dV^ and it satisfies the conservation equation, T^ b — 0, where a | 
denotes covariant differentiation on the hypersurface, which implies the explicit equation 

(au a ) la + Puf a = 0, (42) 

this last being obtained by multiplying the tensor T ab by u a . 

To apply the Hamiltonian formalism we should pass to a spacetime characterization of the energy-momentum 
tensor. In order to do so we define the spacetime energy-momentum tensor as 

T» v = T a V Q l e v b . (43) 

Then, in order to write Tj_j_ and Txi we must first write the energy-momentum tensor in an adapted coordinate system 
{T, X, 6 1 , ■ ■ ■ , 6 d ~ 2 } (see FigureH]), where T and X are the time and radial adapted coordinates, and 9 1 , ■ ■ ■ , 9 d ~ 2 are 
the usual angular coordinates. In more detail, the origin of the new adapted coordinate X is at the shell, and the 
direction of any vector in X alone is the same as that of the outward spacetime normal to the surface, £, and the 
direction of any vector in T alone is the same as that of the tangent vector u. Then, 

= {au^u" + P ih^ + u»u u )} S(X) , (44) 

where W v is the intrinsic metric of the shell written as a bulk d dimensional spacetime tensor, h^ v = — , and 
u M is the velocity of the particles of the shell written as a bulk d dimensional vector. When written as a bulk tensor, 
the energy momentum tensor of the thin shell is proportional to a delta function. Now, T±± is given by T MI/ n M n^, 
which means that T^n^n^ = [au^u^n^riu + P (Wn^riv + u^u v n^n v )] S(X) . Knowing that n^n^ = —1, we have 
that (N- 1 -) 2 — f 2 . So uj_uj_ = "f 2 /f 2 . In the same way, we obtain h^ v n^n v = -1 - (R 2 )/(f 2 ) . We arrive thus at 
Tj__i_ = u (j 2 / f 2 ) S(X) . It is, however, necessary to write the delta function as a function of the radial coordinate r. 
This means we have to return to the original Schwarzschild type coordinates of the metrics (|36|) and (|37| . For that, 
we write dt — u l dT + dr = u r dT : + £ r dX . As we are interested in the energy-momentum tensor on the spacelike 

hypersurface S t , we make dt = and obtain the differential relation dr/dX = £ r — (u r /u t )£ i t = f 2 /"f ■ Taking into 
consideration the relation 6(f(x)) — ^iS(x — Xi)j \f (xi)\ , where are the zeros of a general function f(x), and ' 
denotes differentiation with respect to the argument of the function, we have S(X) = (f 2 /j) S[r — R(i~)] . Plugging 
one equation into the other we get 

T xx = <n&[r ~ R(t)] , (45) 

where a is the energy surface density on the shell, 7 is given in l(40|) . and R(t) is the shell radial function. 

For the other components of the matter tensor, Tu, we have first to define — fe- as the projection vectors 

onto the hypersurface St, where y % are the intrinsic coordinates of St- Then, writing the intrinsic metric of dV^ as 
W = g^ v - ', and knowing that = 0, we arrive at 

T±i = — [(<r + P)(-iV ± /- 2 7 ) Ml + Pm] 5{r - R(t)} , (46) 
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where it, are the components of the velocity vector projected onto £ t . It is then possible to write completely the 
matter constraints as 

H { r ] =^^5[r-R{r)], (47) 

and 




One remark is due, regarding the foliation of spacetime. From the way we have sliced the spacetime, by looking at 
the ADM metric ^ and the fact that we are working with spherical symmetry and static background (cf. |25|) ). we 
see that the shift Lagrange multipliers N l are equal to zero. Therefore, every constraint with components in the £4 
hypersurface Hi will not be relevant to our discussion. In the following we will need neither Tj_j nor TL ™ . 



2. The equations of motion 

By varying the Hamiltonian action, written with the constraints spelt out previously, with respect to g, N- 1 , cp, 
and p r the field equations will be, respectively 



= 0. 



r d 2 dr \ 



F(r) 



P 



2fl 



■T±±, 



d<p 



d-1 



r d - 2 f- 



dr 



+ Np = 0, 



(49) 
(50) 

(51) 
(52) 



where F(r) = (1 — f 2 (r))/r 2 , N = N A -(g rr ) 2 , and p(r) is a redefinition of the canonical momentum radial component 
p r through 



p = r 



„d-2 



£}2 



n 



d-1 



-P, 



(53) 



where £1 is the angular part of the determinant g of the intrinsic metric gtj of the hypersurface S t , with g — g rr r d 2 ft, 
and we have used the fact that the metric is diagonal. In the same way we have redefined the current's time component 



j0 _ r d-2 ^ 2 „-0 



We also have used ([3]) to write 



d-2 



1-/ 



2 \ P 



F 



(54) 



(55) 



One should note that varying with respect to N ± and <p is requiring that the constraints be made equal to zero, 
because both are Lagrange multipliers. That is, in the case of N 1 -, the variation implies H± = 0, and varying with 
respect to ip implies E v = 0. 

Equations (|49| - f52|) are valid for any thin spherically symmetric charged shell in Lovelock gravity coupled to 
Maxwell electromagnetism. One has only to give the thin shell delta-function energy-momentum tensor and integrate 
the equations to find the solutions. Note that the Eq. lj50j) is valid only for a delta function energy-momentum 
tensor. Indeed this equation when integrated gives the equation of motion of the shell itself, but if one has a 
continuous distribution of matter, instead of a thin shell, then within the matter there exists time dependence in the 
metric function, which should be taken into account. In the case of the thin shell we were able to ignore this time 
dependence, since inside and outside the shell it is possible to set up static backgrounds. 
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III. THIN SHELL SOLUTIONS IN A VACUUM BACKGROUND WITH SPHERICAL SYMMETRY 

A. Spherically symmetric vacuum solutions: solutions with black hole character, with naked singularity 

character, and empty solutions 

The vacuum solutions are obtained when we consider the matter action equal to zero, Tj_j_ = 0, and when there is 
no vector current, J M = 0. So, after integrating the field equations above, one obtains the vacuum solutions [13] 

N = N 00 = l, (56) 

2 

f(r) = l+ r ¥ -X9k(r), (57) 

PW = e ^2j' ( 58 ) 
e Q 
(d - 3) r( d " 3 ) 



<p(r) = (/?oo+ , , 6 „, -^r- (59) 



The function gk{r) is defined by 



i/fc 



, { 2G k M + 6 d -2k,i tG k Q 2 
9k(r)={ (d _ 3) r2{d _ k _ 2) j , (60) 

where the integration constants M and Q are the mass and electric charge of the solutions, and where we have used 
Eq. {5} for the definition of k in terms of From Eqs. (|56|) - ([60l) . in particular from Eqs. (|58]) and l|60jl . one notes 
that the vacuum solutions yield a natural division of the Lovelock theory into two branches, namely d — 2k — 1 > and 
d — 2k — 1 = 0, with d being the dimension of the spacetime and k the parameter that gives the highest power of the 
curvature in the Lagrangian. The branch d — 2k — 1 > embodies general relativity when k = 1 (any d), Born-Infeld 
when k = [^y^], and other generic cases. The branch d — 2k — 1 = embodies Chern-Simons type theories alone. 
Worth mentioning now is the fact that in the d — 2k — 1 > branch the empty vacuum solutions have M = and 
Q = 0, whereas in the d — 2k — 1 = branch the empty vacuum solutions have M = — (2Gk)~ 1 and Q = 0. In the 
solutions, there appears an additional unusual parameter \- The parameter \ is the spacetime character, given by 

X=(±l)*+ 1 . (61) 

This means that the vacuum solutions may have a black hole character or a naked singularity character. For x = 1 
the solutions, being of the type found in general relativity, have a black hole character, since for a correct choice of 
the parameters, such as mass and charge, the solution is a black hole solution, although, of course, for other choices 
of parameters it can be an extremal black hole or a naked singularity. For \ = — 1 the solutions, being of a new type 
not found in general relativity, have a naked singularity character, as there is no possible choice of the parameters 
that gives a black hole solution, the full vacuum solution is always singular without horizons. Note also that Eq. 
(|6"T|) implies that if k is odd (such as in general relativity, where k = 1) the character \ = 1 only, whereas if k is 
even (such as in general relativity with a Gauss-Bonnet term, where k = 2) the character x can have both values 
±1. Note also that in the d— 2k — 1 = Chern-Simons theory, the black hole vacuum spacetime, i.e., the vacuum of 
the x = 1 character black hole solution, given by M — and Q — 0, is different from the usual AdS spacetime, and 
there is a mass gap between these spacetimes, with the latter being obtained for M = —(2Gk) , as mentioned above. 
Equations f5"6j) - (|6"Tj) provide the electrically charged solutions of the Lovelock theory defined by the parameters d and 
k coupled to Maxwell electromagnetism, with d > 4. The dimension d = 3 yields singular solutions in the charged 
sector, and so we do not consider this dimension in this work, although in the uncharged sector it yields perfectly 
sensible solutions [l9| . 

Now, given a spacetime solution one should search for horizons and singularities. For x = 1 the solution has a black 
hole character and can have horizons and singularities, i.e., the solutions represent black holes when the parameters are 
appropriately chosen. For ^ = — 1 the solution has a naked singularity character, and should have only singularities 
with no horizons, independently of the choice of the other parameters. First we locate the horizons for the x = 1 
solutions, and then we locate the singularities for both x = ±1- For x = lj the zeros of f(r), in the coordinates used, 
give the horizons. Upon close scrutiny, the properties of these solutions have many similarities with the d dimensional 
Reissner-N6rdstrom-AdS black holes in pure general relativity. Moreover, this set of black hole solutions reduces to the 
d dimensional Reissner-Nordstrom- AdS black holes for k = 1, and the charged Born-Infeld and charged Chern-Simons 
black hole solutions are recovered for d — 2k + 2 and d — 2k + 1, respectively [ijj. For generic values of d and k, in 
analogy with the Reissner-Nordstrom geometry, the black hole solutions possess, in general, two horizons located at 
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the roots of f 2 (r) , one is the event horizon r + , and the other the Cauchy horizon r_ , with r_ < r + . When the two 
horizons merge the black hole is extremal as usual. If the solution is overcharged one has a naked singularity for some 
value of Q in terms of M [53]. Now one finds the location of the singularities, which can be treated together for both 
characters x — ±1- The scalar curvature is 



R = 



1 



2 ^ r 2 



„d-2 



9k{r) 



2 \ 1 



(62) 



For k = 1, general relativity with its Reissner-N6rdstrom-AdS black holes, one finds that R = as it should, thus in 
this case the singularities are located when the Kretschmann scalar, R^ vp<T R^f" , blows up, which is at r = 0. For 



k > 1 the Ricci scalar is no more zero, R ^ 0. In this case, Eq. l(62j) has a singularity at the zero of the function 
<7/c(r), and is due to the existence of an electric field. This is a real timelike singularity located at r e 



2(d-3) (M + (2G fe )-M rf _ 2M ) 



l/(d-3) 



fc > 1 



(63) 



This singularity can be reached in a finite proper time interval. However, an external observer is protected from it 
because both horizons cover it, r e < r_ < r+. Moreover, for even fc, regions where r < r e have a metric with complex 
functions, so there is no solution in this region. For odd fc, the metric can be defined in the region < r < r e , where 
r = is also a spacetime timelike singularity. This solution has no horizons, and is defined between two timelike 
singularities and so is of no interest to us. For all purposes the solutions of interest, for fc > 1, are defined in the 
interval r e < r < oo, and for k = 1 (general relativity) in the interval < r < oo. The corresponding Carter-Penrose 
diagrams can be easily sketched. 



B. Spherically symmetric thin shell solutions through the Hamiltonian formalism 

1. Shell dynamics 



(a) The two master equations 

We now want to find the equations governing the motion of a thin shell in a vacuum background of a given Lovelock 
theory specified by the dimension d of the spacetime and the parameter fc, which gives the highest power in the 
curvature terms in the Lagrangian of the theory. In order to obtain shell solutions in a vacuum background, as 
opposed to pure vacuum solutions, we study the complete field equations, Eqs. (|49|) - ([52"1) . and take into consideration 
the matter term and the electric current J M in those equations. We know that inside the shell the spacetime 

solution, with mass Af_ say, is obtained by integrating the constraint Ti± = 0, i.e., the Eq. lj50j) with T±± = 0, with 
the vacuum solution given in |57l) . In the same way, outside the shell the spacetime vacuum solution, with mass M + 
say, is obtained by integrating the constraint H± = 0, i.e., the Eq. f50|) again with Tj_j_ = 0, with the vacuum solution 
given in f57j) . It remains now to integrate the constraint H± = in the neighborhood of the shell, from R — e to 

R + e, where e is the infinitesimal thickness of the shell, with e — > 0. That is, we have to impose drH±_ = 0. 
This integration should be performed in the asymptotic region, outside of the black hole, since there the normal to 
the t = constant hypersurfaces is timelike, and thus the Hamiltonian formalism is directly applicable. In addition, 
one has to fix from the start the same value of the character \ on both sides of the shell (one has to choose either 
X = 1 or x = ~ l)j i- e -i one has to choose which of the two types of spacetime one is working with. Using then Eq. 
(l50ll we obtain 



R+e 
R-e 



dr —K 



(d-2)\ d 



„d-2 



dr 





i" 











2 



T±x }= , 



(64) 



where T±± is given in (|45|) . Knowing ((3]) and J5]), using {55]), and defining the energy-density content of the shell 
to = a £ld-2R d ~ 2 , the integration yields 



-m( 7+ + 7 _) = (M+-M-) 



<Q\-Q 2 -) 1 
2(d-3) R d ~ 3 ' 



where following Eq. (J3D] 



7+ 



fl + R 2 and 7. 



fi+R 2 . 



(65) 



(66) 
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7± being the generalized Lorentz factor, and where the different indices + and — stem from the fact that the integration 
is made both in the y( - ) (from R — e to R) and V^ + ' (from R to R + e) spacetimes. Note that when integrating 
Eq. (|64| . we have used the fact that the radial electric field is zero inside a uniformly charged sphere, due to Gauss' 
theorem, and that outside it is (Q+ — Q-)/{r 2 + e 2 ), where e denotes the distance from the surface on the outside to 
the point of measurement, r denotes the radius of the sphere, and f~ 2 {r) is continuous in the domain of integration. 
When there are black holes this domain should be exterior to horizon of the black hole solution, where the solutions 

(e) 

are static. Then, the integration in r of the H ± between the limits R — e and R + e is, in the limit of e — > 0, 
equal to zero. Thus, the electric field does not contribute to the shell equation through the radial integration of the 
Hamiltonian constraint. It is sometimes useful to write Eq. lj65|) in another way, namely to multiply both sides by 
7_ — 7+ . This yields the following equivalent equation 

m X (9k+(R)- 9k-(R)) = [9k + (R) h - 9k-(R) k ] i?^ 2 ^ 1 ( 7 - - 7+) ■ (67) 

Eq. (|65l) . and its equivalent (|67|l . were deduced for the asymptotic region, where the vacuum spacetime is static, but 
in the case a black hole is present or being formed, with some care, those equations can be extended, by continuity, 
to the interior black hole region. 

We proceed our study by analyzing now the electrodynamic constraint l|5Tj). Inside the shell the electric field 
solution, Q_/r d_2 , is obtained by integrating the constraint E v = 0, i.e., Eq. lj5l]) with j° = 0, yielding the vacuum 
solution l(58|). In the same way, outside the shell the electric field solution, Q+/r d ~ 2 , is obtained by integrating the 
constraint E v = 0, i.e., Eq. (5lj) again with j = 0, yielding the vacuum solution (|58| . It remains now to integrate 
the constraint E v in the neighborhood of the shell, from R — e to R + e, where e is the infinitesimal thickness of the 

shell, with e — > 0. That is, we have to impose dr E v = 0. Using then Eq. (f53~|l . we obtain 

fR+e j r R+e 

/ dr-f (r d " 2 p)= / drr d - 2 f, (68) 

JR-e dr jR-e 

where j° is the time component of the vector current, and j° = cr e S(r — R)u a , where cr e is the surface charge density. 
Defining the charge of the shell as q = er e O d _ 2 i? <i-2 yields 

Q+-Q-=q. (69) 

Relation 1(69]) is Gauss' law, and is a trivial consequence of the conservation of charge. 

The two master equations are then Eq. (|65|l (or its equivalent l(67j) ) and Eq. |69|). The former can be seen as a 
dynamic equation for R 2 , the latter is a static equation that simply expresses the conservation of charge. Thus in the 
rest of our study we concentrate on Eq. (|65ll (or equivalently on Eq. (|BT|l ). 

(b) Rewriting the equations to simplify the analysis 

Equation f65|) (or Eq. (|67|) ) was derived for the asymptotic flat region, thus if there is a black hole present it was 
derived for the region exterior to the event horizon. We should now transform it in order to have it written in a more 
workable manner. For that we have to square Eq. lj65]l appropriately. To simplify the notation let us define the right 
hand side of Eq. 1(65]) as E, i.e., 

Es( M + -M.)-!<|^f>^. (70) 

Note that from (|65|) one sees that this quantity E is positive for positive shell mass m (i.e, positive energy density 
a > 0), a condition we always assume throughout. Now we obtain from Eq. (|65| a set of two equations which give 
the dynamics of the shell in detail. First, squaring Eq. lj65j) and using the definition lj70j) . we get for R 2 the following 
expression 

2 - f- , (71) 
or, alternatively 

2 

- fl ■ (72) 

Second, there is one last step needed for the entire set of relevant equations, in the asymptotic region, to be written 
down. This has to do with the fact that squaring Eq. I|65p yields Eq. (|7T|l (or, alternatively, (72])) making the 
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solutions of l(65|) also solutions of (|7T|l (or, alternatively, f72|) ) but not all solutions of the (fTTj) (or, alternatively, (|72|) ) 
are solutions of lf65|) . So, a sufficient condition is found by putting back Eqs. ([7Tjl and (|72l) appropriately back into 
Eq. (|65|) . One arrives then at the following conditions 

ui 

fl-f-< -^E 2 . (74) 

Equations f73|) - ([74l) define implicitly a constraint radius r c , meaning that, in the asymptotic region, a solution for 
R(t) is a solution of the equation of motion only if it obeys f7Tj) (or, alternatively, (|72|) ) and R{t) > r c for all R. We 
need now to choose which of the two relations, Eqs. |73|) and |74|l . hold for the system in question. In order to do so, 
we have to know what signal the difference f\ — fl carries. From lj57|) one sees that the most general expression is 
/+ — /_ = X (9k-(r) — cjk+(r)) , where gk(f) comes from ([SO)) . The signal on right hand side of this equation depends 
on the difference x (<7fc-( r ) — 9k+(r)) being positive or negative, and this determines which relation in ([73|) - (|74"1) is the 
relevant one to use. For example, in the case where the interior is flat, <?fc-(r) = 0, the right hand side is ~X9k+{r)- 
In this case when — X3fc+( r ) < the relation in lj73j) is the relevant one, when — xgk+(r) > the relation in (|74l) is the 
appropriate one. (Note that for pure vacuum, no shell, — X9k+{.r) > 0, represents a solution with naked singularity 
character, while — xgk+(r) < represents a solution with black hole character). In brief, Eq. f73|) defines a constraint 
radius r c which is of use in the asymptotically exterior region. So, in the asymptotically exterior region, Eqs. ([Till (or 
([72]l) and (JT3J) - (JT4J) recovers the same solutions as those of the shell equation lj65j) . with lj73|) - (|74|) being a constraint 
equation. With this set of equations, containing no square roots, one can analyze the spacetime evolution of the shell, 
in the asymptotic region, outside the black hole. 

It is also useful to derive the acceleration of the shell radius. Differentiating Eqs. lf7T|l - lf72|) with respect to the 
proper time of the shell r, we get 

mR = m 2 filT m + (d - 2)R d -^Pl + l- 
m 2 ( df 2 df 2 \ 

where we have used the conservation equation lj42j) to arrive at the result ^ = — (d — 2)R d ~ 3 PR, which was used in 



2. Analysis of two cases: matter in equilibrium and gravitational collapse 

We can now apply the above equations of motion to a variety of physical situations. We choose two interesting 
cases, namely static matter in equilibrium and gravitational collapse of matter. 



(a) Matter in equilibrium 

To find a solution for a static shell we have to determine the radii where both the velocity of the shell and its 
acceleration are null, i. e., R = and R = 0, in the Eqs. l(7"lj) (or, alternatively, (|72|) ) and ([75| . respectively. Forcing 
R = in dn]) (or, alternatively, {22])) yields 



R 



d-3 



<Q\ - Q-) 

2(d- 3) 



(76) 



where Rq is the radius at which the shell is static, and /+ and /_ are evaluated at Rq. To be in equilibrium we have 



to impose, in addition, R — 0, which when combined simultaneously with f76j) gives the pressure necessary to hold 
the thin shell in equilibrium. This pressure is then given by 



P 



1 



(d-2)R? ) - i Oi-2 7+7- 



{Q 2 + -Ql)l + l-m 
2ER^ 2 




+ 7+ 



dfl 
dR 



(77) 



where the functions are evaluated at Rq, given by lf76j) . Thus, we obtain the pressure in terms of the parameters of 
the problem, such as d, fc, m, M + , M_, Q + , Q_, and the radius of the static shell Rq, for a static configuration. As 
a particular simple case we may apply the above expressions to find the pressure necessary to hold the shell in static 
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equilibrium in general relativity (k = 1, x = 1)> without charge, and for zero cosmological constant and flat interior. 
This gives 

p (^-3) Gm 2 

2(d-2) n d _ 2 < 2 « 3 -m)' 1 J 

where i?o is the radius of the static shell, given by i?o = (Gm 2 /[2(m — M)]) 1 ^ d 3 \ and G is the Newton's constant 
for k = 1, i.e., we have done Gfc = i = G, whatever the dimension d. In four dimensions, d = 4, Eq. |78|) reduces to 
P = (Gm 2 )/(16 7T^(i? - m)) , with R Q =Gm 2 /[2(m - M)], confirming the result given in (27j . 

(b) Gravitational collapse 

In this analysis of gravitational collapse in Lovelock gravity coupled to Maxwell electromagnetism, we first study 
a simple example of charged dust matter collapsing into an empty interior and then prove cosmic censorship in the 
generic case of shell collapse into an interior free of naked singularities. An analysis of gravitational expansion, not 
worked out here, follows straightforwardly by performing a time reversal operation on the collapsing solutions. 

(i) Gravitational collapse of dust matter in an empty interior: 

In order to have an idea of the main features of gravitational collapse in Lovelock gravity coupled to Maxwell 
electromagnetism, rather than analyzing the full problem, which can be daunting since there are a great number 
of parameters to play with, we investigate the simplest problem, namely, of gravitational collapse of charged dust 
matter into an empty interior. Even then, in this case the task of studying it in full detail is enormous, since one 
has first to set the dimension of the spacetime, the dimension of the k parameter, and then all the other parameters. 
So we resort to study it somewhat generically, and then give some particular examples in the diverse theories that 
we chose to study. The aim is to show that gravitational collapse is possible in this subset of theories of Lovelock 
gravity, reinforcing some similarities which this has with general relativity, as well as showing some new features. So, 
in the following we analyze the case where the matter is composed of dust particles with P — 0, in which case to is a 
constant and can be identified with the shell's rest mass, and where the interior is empty, in which case the function 
f-(R) is given through f 2 (R) = R 2 /l 2 + 1. 

Before proceeding, it is necessary to emphasize that the Lovelock theory we are analyzing is separable into two 
branches: (i) the d — 2k — 1 > 0, which can have even and odd dimensions, and includes general relativity when 
k = 1 (any d), includes the Born-Infeld case of the dimensionally continued theory when k = [^-], and includes other 
generic cases, and (ii) the d — 2k — 1 = 0, which can only have odd dimensions and is precisely the Chern-Simons case 
of the dimensionally continued theory. We study each branch in turn. 

The branch d — 2k — 1 > (general relativity when k = 1 (any d), Born-Infeld when k — [^jp], and other generic 
cases): 

The d — 2k — 1 > branch allows theories in even and odd dimensions. For instance in d — 6 one can have k = 1 
(general relativity) and k = 2 (general relativity with a generalized Gauss-Bonnet term) theories, which in the case 
k = 2 gives a Born-Infeld theory. On the other hand, in d — 7 one can also have theories with k — 1 (general 
relativity) and k = 2 (general relativity with a generalized Gauss-Bonnet term) , but none of these is a Born-Infeld 
theory. Indeed, as already alluded to in Sec. [ill the Born-Infeld case is the realization of this Lovelock theory in even 
dimensions, for which k = 

Considering thus an empty interior, with M_ = and Q- = 0, one has / 2 (i?) = R 2 /l 2 + 1. Putting then M + = M, 
Q+ = Q, and \ = (±l) fe+1 (see JED)) one finds from Eqs. l[71 ^ -([72 | that the shell equation for d - 2k - 1 > is, 

M 



p& = | 2 (rf-3) R*- s , m 



4 M 



2(d-3)fl d " 3 

. 2Gu \ 1/k f n r e Q 2 ^ 1/fcN 2 



Rd-2k-iJ y 2(d-3)R d - 



R \ 

In order to have physical solutions, we have to demand R 2 > 0. 

The equation of the acceleration of the thin shell lf75)l . allows us to understand the forces acting on the shell in 
collapse. Expanding the 7± defined in l|66p . with R 2 replaced by the right hand side of ifflj) or (|72jl . appropriately 
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FIG. 2: In this figure the branch d — 2k — 1 > (general relativity when k = 1, Born-Infeld when fc = [^p], an d other 
generic cases) is studied in some particular instances. In each plot the effective potential R 2 as a function of r (which gives the 
turning points and the allowed regions for the shell), and the metric function f 2 (r) as a function of r (which gives the horizon 
formation, for several different values of the parameters) are displayed. The vertical axis 3>(r) represents both R (r) (full line) 
and f+(r) (dashed line). Note that the thin shell trajectory is allowed only in the region where R 2 is positive. In all the plots 
the respective Newton's constant is put equal to the unity, Gk = 1, so that radii and energies are measured in the respective 
Planck units (h = 1, c = 1). (a) General relativity, k = 1, left plot is for d = 4 and right plot for d = 10, with x = 1 in both 
plots (x = —1 is not possible for k odd), (b) Born-Infeld, k = left plot is for X = 1 an d right plot for x = — 1, (c) General 

relativity with a generalized Gauss-Bonnet term, k = 2, a particular case of 1< k < left plot is for X = 1 and right plot 

for x ~ ~ 1- See text for details. 



chosen, we have 
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Then Eq. l(75|) for an empty interior turns into 

- 2R Q 2 7+ 7- rn_ ( 1 / 2G k M eG k Q 2 \ ^ 

P 2ER d - 2 X 2E X y- k yRd-^- 1 d-iR^ d - k -^)J 

( 2{d-2k-l)G k M 2{d- k-2)eG k Q 2 \ 

where the functions 7± are given in (JHOJ) , and we have used Eq. (|65|) . The first term on the right hand side of Eq. 
(|8l]) is proportional to the radius of the thin shell, meaning that this term dominates the acceleration for large values 
of R, it tends to — oo as R — ► oo dominating all the others, and since the term is negative it points towards the 
center. The term proportional to the square of the electric charge Q 2 includes a R 2 terms and a \R\ term (due to 
the product 7_7 + ). This means that the charge term is positive for a thin shell, works to push the shell to higher 
radii, and it tends to zero as R — > oo. The third term has several terms included, like gravitational and viscosity 
forces per unit mass terms. A simple example where the terms in Eq. f81j) are drastically reduced is when we put 
d = 4, k = 1, I 2 — ► oo, Q = 0, i.e., the case of an uncharged shell collapse in usual general relativity in four dimensions. 
Then we get for the acceleration of the shell, R = —M/R 2 + m 2 /(4i? 3 ) . The term m 2 /(4i? 3 ) is a correction to the 
Newtonian dynamics, representing a self gravitational potential energy. 

Since, in general, the character \ has two values, \ — ±1, for the branch we are analyzing, we study separately both 
cases. Note first that for x — 1 one has from l(57|) that the shell spacetime solution has a black hole character (i.e., 
for certain choices of the parameters one has a black hole solution) , whereas for \ = — 1 one has from ([57)1 that the 
shell spacetime solution has a naked singularity character (i.e., there are no possible choice of parameters that give a 
black hole). So, since when k is odd (such as in general relativity where k = 1) the character \ — 1 for sure, odd k 
allows black hole character solutions only, and since when k is even (such as in general relativity with a Gauss-Bonnet 
term where k = 2) the character x can have both values ±1, even k can have both black hole and naked singularity 
character solutions. 

For x — 1 the spacetime has a black hole character. When the shell collapses it will pass through its horizon 
radius, defined by the equation f+{rh) = 0. So for % = 1 shell collapse implies the formation of a black hole. To 
see this more clearly, note that from Eq. l[57|) the relation f 2 (rh) — defining the horizon turns into (l + r 2 /l 2 ) = 

(2G k /r d i ~ 2k ~ 1 ) ^ (A/ — (eQ 2 )/(2(d~ 3)rf^ 3 )) k . One can show by inspection that the horizon radius is always 
larger than the turning radius r t , r/i > r t . So in this type of Lovelock gravity the bounce is always inside the event 
horizon, which implies the shell expands into another universe, a result already obtained in pure general relativity 

For x = ~ 1 the spacetime has a naked singularity character. We can now give an argument which shows that 
cosmic censorship holds in the case of d — 2k — 1 > 0, with an empty interior. The spacetime would be singular if the 
shell would hit r e , given in {63]). Through the shell equation, Eq. ^B, we see that M - (e Q 2 )/{2 (d - 3)R d ~ 3 ) > . 
However, the singularity at r e makes the right hand side of the shell Eq. f65|) zero. This shows that the radius of this 
singularity is out of the region of validity of the shell equation, preventing the collapse of the charged shell to form a 
naked singularity, and validating here the cosmic censorship hypothesis. Note that this is only applicable to k even. 

In Figure[2]the branch d— 2k — 1 > (general relativity when k = 1, Born-Infeld when k — [=-^1, and other generic 
cases) is studied in some instances. In each plot the effective potential R 2 as a function of R (which gives the turning 
points and the allowed region for the shell motion), and the metric function f 2 (R) as a function of R (which gives 
the horizon formation), for several different values of the parameters, are shown. The vertical axis $(r) in the figure 
represents both R 2 (r) (full line) and f 2 (r) (dashed line). The thin shell trajectory is allowed only in the region where 
R 2 is positive. In all the plots we have made the respective Newton's constant equal to the unity, G k — 1, so that 
radii and energies are measured in the respective Planck units (h = 1, c = 1). We have not plotted R since it yields 
a totally different scale in the vertical axis, rendering the other two functions, R 2 and f 2 (R), almost invisible. We 
now look at each choice of plot in turn, (a) General relativity, k = 1. Left plot: d = 4, I = oo, M = 1000, Q = 300, 
to = 20. Right plot: d = 10, I = oo, M — 200, Q = 300, m = 30 (since k is odd, one has x = 1 for both plots). 
More specifically, the left plot is the usual general relativity with d = 4 (which incidentally is also a Born-Infeld 
type theory), and the right plot is general relativity with d — 10. There is always collapse to a black hole, a result 
confirmed by both plots. Indeed, since R 2 is nonzero at the horizon, where f 2 (r) = 0, the shell passes smoothly 
through the horizon itself. So there is no bounce outside the horizon, the bounce being inside the horizon, and into 
another universe. For example, in d = 4 the value of the radius of the bounce, or turning point, inside the horizon is 
given by r t — (to 2 — Q 2 )/ (2(m — M)) = 45.7 in the units defined above. This agrees with the expressions given in [29| . 
where using Israel's formalism [2l|, as opposed to the Hamiltonian methods used here, a full study of the equation of 
motion and of the turning points in d dimensional general relativity was performed, (b) Born-Infeld, k — [^j^]. Left 
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plot: d = 10, k = 4, x = 1, i = oo, M = 10000, Q = 30000, m = 1. Right plot: d = 10, fc = 4, * = -1, / = oo, 
M = 10000, Q — 30000, m = 1. More specifically, the left plot is for \ = 1, and we find there is always collapse to a 
black hole, since i? 2 is nonzero at the horizon, where f 2 (r) = 0, and the shell passes smoothly through the horizon 
itself. There is a bounce inside the horizon into another universe. The right plot is for \ — ~ 1, with the same choices 
of the other parameters as for % = 1- I n this case the external spacetime has a naked singularity character. Here there 
is no horizon because the metric function f\(r) has no zero. There is a zero of R 2 , which limits the region R 2 > 0, 
the region where the shell has its trajectory. The zero is larger than the radius at which there is a singularity. Thus 
the collapse never forms a naked singularity, (c) General relativity with a generalized Gauss-Bonnet term, k = 2, a 
particular case of 1 < k < [^jp], (here the intermediate cases, (1 < k < [^p]), where k is neither minimum (general 
relativity), nor maximum (Born-Infeld) , are studied). Left plot: d = 10, k = 2, \ = 1, I = °°, M = 1000, Q — 300, 
rn = 20. Right plot: d= 10, k = 2, x = — 1, I = oo, M = 1000, Q = 300, rn = 20. The choice of the coefficients a p are 
given in Eq. J3]). More specifically, the left plot is for X = 1, and we see there is always collapse to a black hole, since 
R 2 is nonzero at the horizon, where f 2 (r) — 0, and the shell crosses the horizon smoothly. There is also a bounce 
inside the horizon. The right plot is for x = — 1> where the external spacetime has a naked singularity character. 
There is no horizon, because there are no zeros of f 2 (r). The region where R 2 > 0, where the shell has its trajectory, 
is limited on the left by the zero of R 2 . The zero of R 2 is larger than the radius at which there is a singularity. There 
is thus no collapse to form a naked singularity. In all the figures (a), (b) and (c) we have shown bouncing solutions 
only, although one could easily find parameters for which the collapsing shell suffers no bounce and goes all the way 
down to the singularity. Note that a general set up for the study of the turning points and the shell's trajectory in 
a Carter-Penrose diagram for the causal structure (with a corresponding careful analysis of the normal to the shell 
along the trajectory), as was done analytically for d dimensional general relativity in [28.[29|. could also be performed 
here case to case, i.e., giving d and k. However, this is beyond the scope of this work. 

The branch d — 2k — 1 = (Chern-Simons): 

The d— 2k — 1 = branch implies d is odd always. Moreover the theory is of Chern-Simons type. The c?— 2k — 1 > 
branch and the Chern-Simons d — 2k — 1 = branch, are quite distinct. First, in this latter branch, in this charged 
setting, there is no general relativity (k — 1), because the theory is defined only for d > 3, and thus k > 2 (however in 
the uncharged case the theory is well defined for d = 3, yield three dimensional general relativity, see [19]). Second, 
for the d — 2k — \ > branch one has that the vacuum mass is M_ = 0, whereas for the d— 2k — 1 = Chern-Simons 
branch the vacuum mass is M_ = — (2Gfe) _1 . From ([571 and (|60|) one finds that the interior for the Chern-Simons 
theory with M_ = -(2G fe )" 1 is characterized by f 2 (R) = R 2 /l 2 + 1 also. 




FIG. 3: In this figure the branch d — 2k — 1 = (Chern-Simons, d always odd) is studied in some particular instances. In each 
plot the effective potential R 2 as a function of r (which gives the turning points and the allowed regions for the shell), and the 
metric function / 2 (r) as a function of r (which gives the horizon formation, for several different values of the parameters) is 
displayed. The vertical axis $(r) represents both R 2 (r) (full line) and /+(r) (dashed line). Note that the thin shell trajectory 
is allowed only in the region where R is positive. In the plots the respective Newton's constant is put equal to the unity, 
Gk = 1, so that radii and energies are measured in the respective Planck units (H = 1, c = 1). Left plot is for k — 2, d = 5 and 
X = L an d right plot is for k = 2, d = 5 and x = — L See text for details 

From Eq. flMJ the shell equation for d - 2k - 1 = 0, M_ = -(2G fc )~\ M+ = M, Q- = 0, Q+ = Q, and 
X= (±l) k+1 is, 



M+(2G fc )- 1 - 2(d _ 6 3 V» 
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(82) 



In order to have physical solutions, we have to demand R 2 > 0. 

Since, in general, \ — ilj we study separately both cases. Again note first that for x = 1 one has from |57j) that 
the shell spacetime solution has a black hole character, whereas for x = — 1 one has from ijBTf that the shell spacetime 
solution has a naked singularity character. So, since when k is even (such as in general relativity with a generalized 
Gauss-Bonnet term, where k — 2 and d — 5) \ can have both values ±1, even k can have both black hole and naked 
singularity character solutions, and since when k is odd (such as k = 3, and d — 7) x — 1 f° r sure, odd k allows 
solutions of black hole character only. 

For x — 1 ( so & can be odd or even) the horizon is located at the radius Th given implicitly by the equation (see 

423)) (! + r l/ 12 ) = ( 2 Gk) 1/k (M+ (2G fe )- 1 - {eQ 2 )/{2{d-i)r d h - z )) 1/k . When the parameters to, M, Q, I yield a 
solution for the collapsing shell then the shell will pass through its own event horizon at r% and form a black hole. 

For x = — 1 (so fc can be even only), the shell solution has a naked singularity character. We can again give an 
argument to show that cosmic censorship in the particular case of Chern-Simons (d—2k—l = 0) with an empty interior 
holds. Indeed, from Eq. J65|) one finds that the following relation holds M + (2 Gfc)" 1 - (e Q 2 )/(2 (d - 3) R d - 3 ) > . 
Noting that r e , given in (|63[). makes the right hand side zero, this means that this singularity is unattainable through 
gravitational collapse. Thus cosmic censorship holds in charged collapse in a Chern-Simons type theory. This is in 
contrast to uncharged collapse in the same theory where naked singularities can form [26| . 

Finally we write the equation of the acceleration of the thin shell, which allows us to understand the forces acting 
on the collapsing shell. Expanding the 7± defined in 1(66]). with R 2 replaced by the right hand side of (|7Tjl or (|72| . 
appropriately chosen, we have 
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As before, the constraint radius marks the lower bound of validity of the shell's equations. Then the Eq. <[75|l for an 
empty interior turns into 
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(84) 



where the functions j± are given in fM)) . and we used Eq. |65l) in the second equality. For uncharged collapse Q = 
the acceleration is always negative, as can be easily checked from Eq. f84|) . In the charged case Q ^ 0, there is 
no definite sign on the right hand side of Eq. l(84|) . it depends on the radius at which the shell is located. Thus, 
re-expansion, after a collapsing phase, is allowed. 

In Figure [3] we show the effective potential R 2 as a function of R (which gives the turning points and the allowed 
regions for the shell), and the metric function f 2 {R) as a function of R (which gives the horizon formation), for several 
different values of the parameters. The vertical axis <f>(r) in the figure represents both R 2 (r) (full line) and f\{r) 
(dashed line). The thin shell trajectory is allowed only in the region where R 2 is positive. In all the plots we have 
made the respective Newton's constant equal to the unity, G k — 1, so that radii, energies and forces are measured in 
the respective Planck units (h = 1, c = 1). Left plot: d = 5, k = 2, % = 1, I = 50, M = 20, Q = 610, to = 5. Right 
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plot: d = 5, k = 2, x = — 1, / = 50, M = 20, Q = 300, m = 5. More specifically, the left plot shows a region limited 
by two zeros of R 2 , the region where a trajectory of a thin shell is possible. The horizon, given by the larger zero 
of f 2 {R), is inside the region where R 2 > 0, which means that the shell passes smoothly through the horizon. The 
shell then suffers a bounce inside the horizon into another universe. The plot on the right shows a naked singularity 
character exterior spacetime, which means that there is no zero of f 2 (R), hence no horizon. There is, however, no 
collapse to a naked singularity, because the region where R 2 > is limited by two zeros, the smaller being larger than 
the constraint radius. Thus the shell bounces between two extreme values of R, and so does not collapse at all. In 
the figure we have shown bouncing solutions only, although as before one could produce totally collapsing solutions. 
As in the d — 2k — 1 > branch, a general set up for the study of the turning points and the shell's trajectory in a 
Carter-Penrose diagram for the causal structure could be performed here, but this will not be done. 

(ii) Cosmic censorship: 

Given the experience we have acquired with the above examples, one can now study cosmic censorship directly from 
Eq. ([651 or (|67l) . Cosmic censorship holds if no naked singularity forms from gravitational collapse in an initially 
nonsingular spacetime or if no naked singularity forms from gravitational collapse in a spacetime initially containing 
a black hole. We assume Q+ ^ and Q- ^ 0. There are three such cases in our study: first, both spacetime regions 
on each side of the collapsing shell have character x = L an d the interior solution is a black hole spacetime; second, 
both spacetime regions on each side of the collapsing shell have character x — 1 , an d the interior solution is an empty 
vacuum spacetime; and third, both spacetime regions on each side of the collapsing shell have character x — ~ 1 5 an d 
the interior solution is an empty vacuum spacetime. On all these cases cosmic censorship holds. For the first case, 
that both spacetime regions on each side of the collapsing shell have character x = 1 , and the interior solution is a 
black hole spacetime, given the above dynamic equation one can directly prove that the collapse of a charged shell in 
such a background never yields a naked singularity spacetime. To start with we have a shell collapsing into an interior 
black hole spacetime. This means that the x term on the left of l|67l) is positive (otherwise the interior would not be 
that of a black hole spacetime). From Eq. (|65"1) it is known that gk+(R) k — cjk-{R) k > 0, because this expression is 
the right hand side of (|65l) . If then the shell is collapsing onto an existing black hole, and it is such that the result 
would be a naked singularity, i.e., the collapsing shell is sufficiently overcharged, then f 2 (rh) = for a certain r^, 
and /+ (r) would always be larger than zero (because in the exterior spacetime there would be no horizon) . It is then 
clear that the term (7- — 7+) is negative if the shell reaches r^. However, as we are working in the asymptotically 
outside region, then gk{f) > 0. This means that the signal of <?fc(r) is the same as that of gk{r) k , and so the signal 
of gk+{R) k — gk-(R) k is the same as that of gk+(R) — gk-{R)- This results in the fact that the left hand side of 
([67|l is positive at the horizon of the inner black hole, defined by fl(rh) = 0, and that the right hand side of the 
same Eq. (|67|) is negative, which implies a contradiction. For the second case, that both spacetime regions on each 
side of the collapsing shell have character x = 1, and the interior solution is an empty vacuum spacetime, one can 
treat it as a limiting case of the first case, as the mass of the interior black hole goes to zero, or directly using similar 
arguments as above, with the result that a collapsing shell, if undercharged forms a black hole, if overcharged has 
a bounce back from where it came. For the third case, that both spacetime regions on each side of the collapsing 
shell have character x — ~ lj an d the interior solution is an empty vacuum spacetime, one can write Eq. (|67| as 
m = —gk+(R) h ^ 1 i? d ^ 2fc ^ 1 (j_ —7+) . One sees that when the shell approaches r e , its right hand side approaches zero, 
whereas the left hand side approaches a finite value. So there is a bounce and the shell never collapses to a singularity. 
These three cases prove the result, and yield cosmic censorship in charged backgrounds in the Lovelock theory we 
are studying coupled to Maxwell electromagnetism. We now briefly comment on the uncharged case Q+ = Q- = 0. 
One can prove that in this particular case, in the d — 2k — 1 > branch one has no formation of naked singularities, 
whereas in the d — 2k — 1 = uncharged branch there is formation of naked singularities, so that one can say that 
electrical charge acts really as a cosmic censor. 

Thus, a charged shell never develops a naked singularity. This is rather like having a shell with some angular 
momentum, a situation which is much harder to study directly in theories with d > 3. Nonetheless, the special case 
of an uncharged collapsing shell in d = 3 with angular momentum was studied in [2~i.[25|. with the result that indeed 
the angular momentum also prevents the formation of a naked singularity. 

IV. CONCLUSIONS AND PHYSICAL IMPLICATIONS 

There are two main conclusions from this work. One conclusion is that the Hamiltonian formalism is a powerful 
method to treat in a unified way spacetimes composed of several pieces, such as several vacua and thin shells, in 
theories much more complicated than general relativity such as the subset of theories derived from Lovelock gravity 
coupled to Maxwell electromagnetism we have studied. The other conclusion is that, when the spacetimes in question 
have the same character of those spacetimes provided by general relativity (x = 1), the collapse of the thin shells in 
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the backgrounds, black hole or otherwise, of each different type of Lovelock theory is in many ways similar to the 
collapse in general relativity itself, and when the spacetimes in question have the opposite character (x = — 1), some 
other new features appear. This in turn has the following physical implications: if indeed there are extra dimensions 
with a relative large size, as exposed in the Introduction, then the new and the old features, when confronted with 
experimental data, can provide the signature to the uncovering not only of the actual spacetime dimension d, but also 
of the value of the parameter k, i.e., of which particular Lovelock gravity nature picks up at the appropriate scales, 
whether be it general relativity (k = 1, any d), Born-Infeld (k — maximum, even d), Chern-Simons (k — maximum, 
odd d), or other generic gravity (other k, any d). Of course, a full quantum treatment, or even a semiclassical 
approximation, would be much more appropriate for this kind of questions, but for Lovelock type gravities the 
technical difficulties are exponentiated easily, and thus it is advisable to start up with a classical analysis, as we did 
here. 
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